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In addition to its importance in describing high-energy processes themselves, the 
dynamics of multiparticle production is part of the general field of non-linear phe- 
nomena and complex systems. Multiparticle dynamics is one of the rare fields of 
physics where higher-order correlations are directly accessible in their full multi- 
dimensional characteristics under well-controlled experimental conditions. Mul- 
tiparticle dynamics, therefore, is an ideal testing ground for the development of 
advanced statistical methods. 

Higher-order correlations have, indeed, been observed as particle-density fluctua- 
tions. Approximate scaling with finer resolution provides evidence for a self-similar 
correlation effect. Quantum-Chromodynamics branching is a good candidate for 
a dynamical explanation of these correlations in e+e — collisions at CERN/LEP 
and, as expected, also of those in pp collisions at future CERN/LHC energies. 
However, other sources such as identical-particle Bose-Einstein interference effects 
also contribute. 

A particular question at the moment is the smooth transition from the QCD 
branching domain (gluon interference before hadronization) to the Bose-Einstein 
domain (identical-pion interference after hadronization). Both mechanisms have 
clearly been observed in e + e~ collisions at CERN/LEP energies. The large amount 
of high-resolution data being collected at LEP will allow the study of the genuine 
(i.e. non-trivial) higher-order correlations in both domains. 



1 Introduction 

Recent years have witnessed a remarkable amount of experimental and theoret- 
ical activity in search of scale-invariance and fractality in multihadron produc- 
tion processes. In addition to being an important part of high energy physics 
itself, the dynamics of multi-particle production in collisions of elementary par- 
ticles at high energies (multiparticle dynamics) is part of the general field of 
non-linear phenomena and complex systems. Studies of classical and quantum 
chaos, non-equilibrium dissipative processes, random media, growth phenom- 
ena and many more processes have all contributed in revealing the pervasive 
importance of self-similarity, power-laws, and fractals in nature. Research in 
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these fields is still in full swing and continues to uncover intriguingly simple 
and often surprisingly universal behavior in complex, non-linear systems. 

While considerable experience already exists in many fields for the study of 
two-component correlation, it is often in higher-order (i.e. multi-component) 
correlations that the most interesting properties manifest themselves, in the si- 
multaneous interplay of a large number of components. The special significance 
of multiparticle dynamics for the development of advanced statistical methods 
lies in the fact that it is one of the rare fields of physics where higher-order cor- 
relations are directly accessible in their full multi-dimensional characteristics 
under well controlled experimental conditions. 

Higher-order correlations have recently been observed as particle-density 
fluctuations in cosmic ray, nucleus-nucleus, hadron-hadron, e + e~ and lepton- 
hadron experiments. To study these fluctuations in detail, normalized corre- 
lation integrals are being analyzed in phase-space domains of ever decreasing 
size. Approximate scaling with decreasing domain size is now observed in all 
types of collision, giving evidence for a correlation effect self-similar over a 
large range of the resolution (called intermittency, in analogy to a statistically 
similar problem in spatio-temporal turbulence). 

Parton branching of Quantum-Chromodynamics predicts the type of cor- 
relations observed in e + e _ collisions at CERN/LEP and, as expected, also in 
pp collisions at CERN/LHC energies. However, other sources such as Bose- 
Einstein interference of identical particles also contribute. Fast development 
of the applied technology has taken place over the last few years, in particular 
in the extension of originally one-dimensional to full three-dimensional phase 
space analysis. 

2 Methodology 

2.1 Particle Densities 

A collision between particles a and b is assumed to yield exactly n particles 
in a sub- volume f2 of total phase space O a n . The single symbol y represents 
the kinematical variables needed to specify the position of each particle in this 
space (for example, y can be the full four-momentum of a particle and fl a 
cell in invariant phase space or simply the cm. rapidity]^ of a particle and f2 
an interval of length 5y). The distribution of points in SI can then be char- 
acterized by continuous probability densities P n {yi, ■ ■ • ,Vn)\ n = 1,2, . . .. For 
simplicity, we assume all final-state particles to be of the same type. In this 

6 Rapidity y is defined as y = ^ ln[(E + pi J )/(E — Pl)], with E the energy and pl the 
longitudinal component of momentum vector p along a given direction (beam-particles, jet- 
axis, etc.); pseudo-rapidity is defined as rj = ^ ln[(p + Pl)/(p — Pl)]- 
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case, the exclusive distributions P n (yi, ■ ■ ■ , Un) can be taken as fully symmet- 
ric in yi, . . . ,y n ; they describe the distribution in f2 when the multiplicity is 
exactly n. 

The corresponding inclusive distributions are given for n = 1,2,... by: 
PniVl, ■ ■ ■ ,Vn) = Pn(yi,---,Vn) 

+ Y] — r / P n +m(yu ■ ■ ■ ,y n ,y[, ■ ■ ■ ,y' m ) TT d ^ • ( x ) 

^— m Jo J - J - 

m=l 1/52 i=l 

The inverse formula is 

P n (yi,...,y n ) = p„(yi,...,y n ) 

Here, p n (yi, ■ ■ ■ ,y n ) is the probability density for n points to be at y\,...,y n , 
irrespective of the presence and location of any further points. Integration over 
an interval Q in y yields 

Pi(y)dy = (n) , 

P2(yi,y2)dyidy 2 = (n(n - 1)) , 



a Jn 



dyi... dy q p q (yi,...,y q ) = (n(n-l) ...(n- q + 1)) , (3) 
o Jn 

where the angular brackets imply the average over the event ensemble. 



2.2 Cumulant Correlation Functions 

Besides the interparticle correlations we are looking for, the inclusive g-particle 
densities p q (jji, ■ ■ ■ , y q ) in general contain "trivial" contributions from lower- 
order densities. It is, therefore, advantageous to consider a new sequence of 
functions Cq(y\, . . . , y q ) as those statistical quantities which vanish whenever 
one of their arguments becomes statistically independent of the others. Devi- 
ations of these functions from zero shall be addressed as genuine correlations. 

The quantities with the desired properties are the correlation functions 
- also called (factorial) cumulant functions - or, in integrated form, TJaiele's 
semi-invariantsBI A formal proof of this property was given by KuboB The 
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cumulant correlation functions are denned as.in-.the cluster expansion familiar 
from statistical mechanics via the sequence: IjoO 

Ml) = C X {1), (4) 
pa(l,2) = C 1 (1)C 1 (2) + C 2 (1,2), (5) 
p 3 (l,2,3) = Ci(l)C 1 (2)Ci(3) + Ci(l)C 2 (2,3) + Ci(2)C 2 (l,3) + 

+ C 1 (3)C 2 (1 ) 2) + C 3 (1,2,3); (6) 

and, in general, by 

ftn(l n>) = E E [Ci()---C 1 ()][C 2 (,)...C 2 (,)].-. 

m ii factors i 2 factors 

■■■[C m {,...,)---C m {,. ..,)]. (7) 
v „ ' 

l m factors 

Here, k is either zero or a positive integer and the sets of integers {k} m satisfy 
the condition 



ili — m. (8) 



The arguments in the Ci functions are to be filled by the m possible momenta 
in any order. The sum over permutations is a sum over all distinct ways of 
filling these arguments. For any given factor product there are preciselyQ 

(9) 



[(l!)ii(2!)fe---(m!HyW"-W 
terms. 

The relations (^) may be inverted with the result: 

C 2 (l,2) = p2(l,2)-pi(l)pi(2), 
C 3 (l, 2, 3) = P3 (l, 2, 3) - W2, 3) + 2pi(l)pi(2)p x (3) , 

(3) 

C 4 (l,2,3,4) - /94 (l,2 ) 3,4)-^p 1 (l)p 3 (l ) 2,3)-^p 2 (l,2)p 2 (3,4) 

(4) (3) 

+ 2 £> x (1)^(2)02(3, 4)- 6^(l)pi(2)pi(3)pi(4). (10) 

(6) 

In the above relations we have abbreviated C q {y\, . . . , y 9 ) to C g (l, 2, . . . , g); 
the summations indicate that all possible permutations must be taken (the 



number under the summation sign indicates the number of terms). Expressions 
for higher orders can be derived from the related formulae given inu. 

It is often convenient to divide the functions p q and C q by the product of 
one-particle densities, which leads to the definition of the normalized inclusive 
densities and correlations: 

r q {yi,---,y q ) = p q (yi,---,y q )/pi(yi).-.pi{y q ), (11) 
K q (yi,...,y q ) = C q (y y q )/p 1 (y 1 )...p 1 {y q ). (12) 

In terms of these functions, correlations have been studied extensively for q = 2. 
Results also exist for q = 3, but usually the statistics (i.e. number of events 
available for analysis) are too small to isolate genuine correlations. To be able 
to do that for q > 3, one must apply moments defined via the integrals in 
Eq. 3, but in limited phase-space cells. 



2.3 Cell- Averaged Factorial Moments and Cumulants 

In practical work, with limited statistics, it is almost always necessary to per- 
form averages over more than a single phase-space cell. Let VL m be such a cell 
(e.g. a single rapidity interval of size 5y) and divide the phase-space volume 
into M non-overlapping cells tt m of size 5£l independent of m. Let n rn be 
the number of particles in cell Q m . Different cell-averaged moments may be 
considered, depending on the type of averaging. n 

Normalized cell-averaged factorial moments El are defined as 

771 — 1 V ' 



j_ A h y Pg(yiT-->yg)ULi d yi 
mJL (h y p(y)dy) q 



m— 1 

M 



(14) 



\p / Pg(y\,---,yq)T[ q i=idyi n ^ 

M{5yY^J 5y {p m f 



The full rapidity interval AY is divided into M equal bins: AY" = M 5y; each 
yi is within the <5y-range and (n m ) = ~p m Sy = J s p\{y)dy. An example for 
q = 2 is given in Fig. la. 

Likewise, cell-averaged normalized factorial cumulant moments may be 
defined as 

k <R«\ 1 f c q(yi>--->yg)Ui=i d yi n ^ 

Sy 
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Figure 1: Integration domains for a) the second-order factorial moment and b) the second- 
order correlation integral. The asterisk in a) indicates the position of a particle pair with 
1 3/1 — 3/2 < &y that is excluded from the F2 calculation due to the binning. In b) the asterisk 
is included. 



They are related to the factorial moments by [] 

F 2 = l+K 2 , 

F 3 = 1 + 3K 2 + Kj_ , 

Fi = 1 + 6X2 + 3X1+4^3 + ^4- (17) 

In i*4 and higher-order moments, "bar averages" appear defined as AB = 
J2A m B m /M. 

m 

To detect dynamical fluctuations in the density of particles produced in 
a high-energy collision, a way must be devised to eliminate, or to reduce as 
much as possible, the statistical fluctuations (noise) due to the finiteness of the 
number of particles in the counting cell(s). This requirement can to a large 
extent be satisfied by studying factorial moments and forms the basis of the 
factorial moment technique, known in optics, but rediscovered for multi-hadron 
physics inQ. This crucial property does not apply to, e.g., ordinary moments 
(n q )/(n) q . 

The property of Poisson-noise suppression has made measurement of fac- 
torial moments a standard technique, e.g. in quantum optics, to study the 
statistical properties of arbitrary electromagnetic fields from photon-counting 
distributions. Their utility was first explicitly recognized, for the single time- 
interval case, km and later generalized to the multivariate case inE3. 

c The higher-order relations can be found ini. 
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2.4 Density and Correlation Integrals 

A fruitful recent development in the study of density fluctuations is the correla- 
tion strip- integral methodEil By means of integrals of the inclusive density over 
a strip domain of Fig. lb, rather than a sum of the box domains of Fig. la, one 
not only avoids unwanted side-effects such as splitting of density spikes, but 
also drastically increases the integration volume (and therefore the statistical 
significance) at a given resolution. 

Let us consider first the factorial moments F q defined according to Eq. 14. 
As shown in Fig. la for q — 2, the integration domain Qb — Sm=i ^ m consists 
of M g-dimensional boxes £l m of edge length Sy. A point in the m-th box 
corresponds to a pair (2/1,2/2) of distance \yi — y 2 \ < Sy and both particles in 
the same bin m. Points with \yi — 2/2 1 < Sy which happen not to lie in the same, 
but in adjacent, bins (e.g., the asterisk in Fig. la) are left out. The statistics 
can be approximately doubled by a change of the integration volume £Ib to 
the strip-domain of Fig. lb. For q > 2, the increase of integration volume 
(and reduction of squared statistical error) is in fact roughly proportional to 
the order of the correlation. The gain is even larger when working in two- or 
three-dimensional phase-space variables. 

In terms of the strips (or hyper-tubes for q > 2), the density integrals 
become 

F q (<%) = 7 — —j—^ / vFf A • (18) 

These integrals can be evaluated directly from the data after selection of a 
proper distance measure — yj\, [{y% — yj) 2 + {<fti — (j^j) 2 ] 1 ^ 2 , or better yet, 
the four- momentum difference Q 2 j — —{vi ~^Pj) 2 ) and after definition of a 
proper multiparticle topology (GHP integrally snake integrally star integral 
E3). Similarly, correlation integrals can be defined by replacing the density p q 
in Eq. 18 by the correlation function C q . 

The numerator of the factorial moments F q can be determined by counting, 
for each event, the number of g-tuples that have a pairwise Q 2 ^ smaller than 
a given value Q 2 and then averaging over all events. Using the Heaviside 
unit-step function 0, this can be mathematically expressed as 

W)=^6 ! e n viQ'-Qk.j) . as) 



il<...<i q aU pairs 



where the factor ql takes into account the number of permutations within a 
q-tuple. 
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The normalization is obtained from "mixed" events constructed by random 
selection of tracks from different events in a track pool. The multiplicity of 
a mixed event is taken to be a Poissonian random variable, thereby ensuring 
that no extra correlations are introduced. A correction factor is applied for 
the difference in average multiplicity of the Poissonian and the experimental 
distribution. The mixed events are treated in the same way as real events. 

2.5 Power-Law Scaling 

The technique proposed ini consists of measuring the dependence of the nor- 
malized factorial moments (or correlation integrals) F q (Sy) as a function of the 
resolution Sy. For definiteness, 5y is supposed to be an interval in rapidity, 
but the method generalizes to arbitrary phase-space dimension, as occurs with 
the use of Q 2 . 

As pointed out above, the scaled factorial moments enjoy the property of 
"noise-suppression". High-order moments further act as a filter and resolve 
the large n m tail of the multiplicity distribution. They are thus particularly 
sensitive to large density fluctuations at the various scales Sy used in the anal- 
ysis. n 

As proven inQ, a "smooth" (rapidity) distribution, which does not show 
any fluctuations except for the statistical ones, has the property of F q {8y) being 
independent of the resolution Sy in the limit Sy — > 0. On the other hand, if 
dynamical fluctuations exist and P p is "intermittent", the F q obey the power 
law 

F q {Sy) tx (8y)-+< , {Sy -> 0). (20) 

The powers <p q (slopes in a double-log plot) are relatedQ to the anomalous 
dimensions 

d q = cf >q /(q-l) , (21) 

a measure for the deviation from an integer dimension. Equation 20 is a scaling 
law since the ratio of the factorial moments at resolutions L and I 

R = F q {t)/F q {L) = {L/C)+* (22) 

only depends on the ratio L/t, but not on L and t, themselves. 

As pointed out above, the experimental study of correlations is already 
difficult for three particles. The close connection between correlations and 
factorial moments offers the possibility of measuring higher-order correlations 
with the factorial moment technique at smaller distances than was previously 
feasible. Via Eq. 21, the method further relates possible scaling behavior of 
such correlations to the physics of fractal objects. 
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Figure 2: Correlation-integral method applied to NA22 hadron-hadron data in terms 

of a) In Fg and b) \nK* as a function of — InQ 2 . 

One further has to stress the advantages of factorial cumulants compared 
to factorial moments, since the former measure genuine correlation patterns, 
whereas the latter contain additional large combinatorial terms which mask 
the underlying dynamical correlations. 

The definition of "intermittency" given in (20), has its origin in other 
disciplines.^] It rests on a loose parallel between the high non-uniformity of the 
distribution of energy dissipation, for example, in turbulent intermittency and 
the occurrence of large "spikes" in hadronic multiparticlc final states. In the 
following we use the term "intermittency" in a weaker sense, meaning the rise 
of factorial moments with increasing resolution, not necessarily according to a 
strict power law. 



3 The State of the Art 

The suggestion that normalized factorial moments of particle distributions 
might show power-law behavior has spurred a vigorous experimental search 

d For a masterly expose of this subject seeH. 
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Figure 3: Comparison of density integrals for q = 2 in their differential form AF2 (in 
intervals Q 2 ,Q 2 +dQ 2 ) as a fup«tion of 2log(l/Q 2 ) for e+e — (DELPHI, preliminary) and 
hadron-hadron collisions (UA1)E3 

for (more or less) linear dependence of \nF q on — ln<5y. A full review of the 
present situation is given inlla. 

As an example, we give Fig. 2a where hadron-hadron data0 on F q are 
plotted as a function of — In Q 2 , with all two-particle combinations in an n- 
tuple having Q 2 j < Q 2 . The following observations can be made: 

i) the moments show a steep rise with decreasing Q 2 ; 

ii) negatives are much steeper than all-charged, 

iii) F| is flatter for (H — ) than for all-charged or like-charged combinations. 
The last two observations directly demonstrate the large influence of iden- 
tical particle correlations on the factorial moments and their scaling behavior. 
These results agree very well with—results from the UA1 collaborationtj In 
Fig. 2b it is, furthermore, shown lia in terms of the K*(Q 2 ) (using the star 
topology) that genuine correlations exist and increase with improving resolu- 
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tion at least up to order q = 5 in hadron-hadron collisions. Again, like-charged 
particle combinations increase faster than all-charged ones. 

Of particular interest is a comparison of hadron-hadron to e + e _ results in 
terms of same and opposite charges of the particles involved. This comparison 
has been done for UA1 and DELPHI data incj and is shown in Fig. 3 for 
q = 2 (in fact, in this figure a differential form AF2 of Eq. 19 is presented). 
An important difference between UAl and DELPHI can be observed in a 
comparison of the two sub-figures: 

For relatively large Q 2 (> 0.03 GeV 2 ), where Bose-Einstein effects do not 
play a major role, the e + e~ data increase much faster with increasing —2 logQ 2 
than the hadron-hadron results. For e + e~, the increase in this Q 2 region is 
very similar for same and for opposite-sign charges. 

At small Q 2 , however, the e + e _ results approach the hadron-hadron re- 
sults. 

The authors conclude that for e + c~ annihilation at LEP at least two pro- 
cesses are responsible for the power-law behavior: Bose-Einstein correlation 
following the evolution of jets. In hadron-hadron collisions at present collider 
energies, the Bose-Einstein effects seem mostly relevant. 

The— exact functional form of AF2 or AK2 is derived from the data of 
UAl Bt3 and NA22P again in its differential form, in Fig. 4. Clearly, the 
low-Q data favor a power law in Q over a (non-scaling) exponential, double- 
exponential or Gaussian law. 

In Fig. 5, the NA22 results for two-, three- and four-particle Bose-Einstein 
correlations^equivalent to AF^Q 2 ) are compared to a multi-Gaussian para- 
metrizationEJ In the conventional linear presentation of Fig. 5a the fits look 
reasonable. If the same data and same fits are presented in a log-log scale, 
however, the power-law deviation from the multi-Gaussians starts to become 
visible (Fig. 5b). 

If the observed effect is real, it supports a view recently developed incll. In 
that paper, intermittency is explained from Bose-Einstein correlations between 
(like-sign) pions. As such, Bose-Einstein correlations from a static source are 
not power-behaved. A power law is obtained if i) the size of the interaction 
region is allowed to fluctuate, and/or ii) the interaction region itself is assumed 
to be a self-similar object extending over a large volume. Condition ii) would 
be realized if parton avalanches were to arrange themselves into self-organized 
critical statesE3 Though quite speculative at this moment, it is an interesting 
new idea with possibly far-reaching implications. We should mention also that 
in such a scheme intermittency is viewed as a final-state interaction effect and 
is therefore not troubled by hadronization effects. 

In perturbative QCD, on the other hand, the fractal structure of jets fol- 
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Figure 4: Density integral AF2 as a function of log Q 2 for NA22EJH-icorrelation integral A.K2 
as a function of log Q for like-charged pairs in UA1 (preliminary) compared to power-law, 
exponential, double-exponential and Gaussian fits, as indicated. 




0.5 1 1.5 2 4 6 

Q 2 (GeV 2 ) -In Q 2 (GeV 2 ) 



Figure 5: The normalized two-, three- and four-pa*Mjicle inclusive densities (equivalent to 
AF g ) as a fimction of Q 2 (left) and — InQ 2 (right U 2 1 Curves show the multi-Gaussian fits 
according toc3. 
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lows, in principle, from parton branching and the intermittency indices <p q 
tly related to the anomalous multiplicity dimension 70 = (Gcts/ft) 1 ' 2 
and, therefore, to the running coupling constant a s via the simple 



relation, 



2 1 

(q - l)D - 2 7o 



where D is the usual topological dimension .of the analysis. In the same the- 
oretical context, it has been argued E3c3 ! Ej that the opening angle between 
particles is a suitable and sensitive variable to analyze and well suited for 
these first analytical QCD calculations of higher-order correlations. It is, of 
course, closely related to Q 2 . 

A first analytical QCD calculation is based on the double-log approx- 



imation (DLA) withpangular ordering EJ (for a recent experimental study of 
angular ordering seeEil) and op_Jocal parton- hadron dualityE3 A preliminary 
comparison with DELPHI dataEj gives encouraging results, but shows that the 
situation is far from trivial, since finite-energy effects, four-momentum cut-offs, 
resonance decays etc. still dominate at LEP energies. 



4 Summary 

Multiparticle production in high-energy collisions is an ideal field to study 
genuine higher-order correlations. Methods also used in other fields are being 
tested and extended here for general application. Indications for genuine, ap- 
proximately self-similar higher-order correlations are indeed found in hadron- 
hadron collisions, but still need to be establised in their genuine and self-similar 
character in e + e~ collisions at high energies. At large four-momentum distance 
Q 2 , they are not only expected to be an inherent property of perturbative QCD, 
but are directly related to the anomalous multiplicity dimension and, therefore, 
to the running coupling constant a s . At small Q 2 , the QCD effects are comple- 
mented by Bose-Einstein interference of identical mesons carrying information 
on the unknown space-time development of particle production during the col- 
lision. The interplay between these two mechanisms, particularly important 
for an understanding of the process of hadronization, is completely unknown 
at the moment. 
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